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^ ABSTRACT 

O ; We present H l synthesis observation of the nearby dwarf galaxy 

to IC 2574 (a member of the M81 group of galaxies) made with the NRAO 

Very Large Array in its B-, C- and D-array configurations at high 



spatial and velocity resolution (95 pc x 2.6 kms ^). In addition, we 



^ \ present optical broad and narrow band images obtained with the 2.2- 

m telescope of the Calar Alto observatory. The VLA H I observations 
show a stunning amount of detail in the form of 48 mostly expanding H I 
shells and holes in its neutral interstellar medium. These features range 
in size from about 100 (a limit set by the size of the beam) to about 1000 
pc dominating the appearance of the H I surface brightness map. Their 
dynamics clearly influence the velocity field of IC 2574. In addition to 
well defined holes, some large scale (> 1000 pc) coherent features are 
visible in the channel maps. They are probably the remainder of an 
older shell population. 



^Visiting Astronomer, German-Spanish Astronomical Centre, Calar Alto, operated by the Max-Planck-Institute 
for Astronomy, Heidelberg, jointly with the Spanish Commission for Astronomy. 
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Current star formation, as traced by Ha emission, is predominantly 
found along the rims of the larger H I holes, suggesting propagating star 
formation. On linear scales of 95 pc, star formation occurs if the 
H I surface density reaches values higher than 10^^ cm~^. The radial 
expansion of the H i holes (8-12 kms~^), the indicative ages (10-60 
10^ yr) and the energy requirements for their formation (10^^-10^^ ergs) 
can be understood in terms of the combined effects of stellar winds and 
multiple supernova explosions of the most massive stars formed during 
a recent phase of active star formation. 

The scaleheight of the H i layer is found to be ^ 350 pc, consider- 
ably thicker than that in massive disk galaxies. This is due to a lower 
gravitational potential (for the same observed onedimensional velocity 
dispersion of ^ 7 kms~^). This puffed up disk implies a lower H i vol- 
ume density (0.15 cm~^). This, combined with the reduced gravitational 
pull and solid body rotation throughout the disk explains why the diam- 
eter distribution in dwarf galaxies extends to substantially larger values 
than in spiral galaxies. 

A comparison with other galaxies shows that the energies needed to 
create these structures are the same for all types of galaxies, at least to 
first order. The overall statistical properties of the H i holes and shells 
in galaxies show clear trends with Hubble type (or rather mass), such 
as in their diameter distribution, expansion velocities and ages. 

Subject headings: Galaxies: individual (IC 2574) - Galaxies: ISM - ISM: 
structure - Radio lines: galaxies - Radio lines: ISM 



2 



1. Introduction 

Some twenty years ago, Heiles (1979, 1984, 
1990) discovered large, shell-like structures 
in the neutral hydrogen (H I) distribution 
of our Galaxy by inspection of maps made 
in the 21-cm hne. Since this discovery, a 
wealth of observations obtained with power- 
ful Synthesis Radio Telescopes such as the 
Very Large Array (VLA) or the Westerbork 
Synthesis Radio Telescope (WSRT) revealed 
similar structures in our nearest neighbors. 
Prominent examples are M 31 (Brinks 1981, 
Brinks & Bajaja 1986), M 33 (Deul & den 
Hartog 1990), Holmberg II (Puche et al. 1992, 
hereafter referred to as PWBR92), the galax- 
ies M 101 and NGC 6946 (Kamphuis 1993) 
and IC 10 (Wilcots & Miller 1998). The lat- 
est, impressive additions to this list are the 
SMC (Staveley-Smith et al. 1997) and LMC 
(Kim et al. 1998). All these observations in- 
dicate that the interstellar medium (ISM) of 
medium to late-type galaxies is dominated 
by features, which are variously described as 
shells, rings, holes, loops, bubbles or cavities. 

Obviously, one of the basic questions is 
which physical process is able to produce 
these observed structures. The general pic- 
ture which has emerged is that they are the re- 
sult of combined stellar winds and supernova 
explosions produced by young stellar associa- 
tions. For review articles on that topic see 
Tenorio-Tagle & Bodenheimer (1988), and 
van der Hulst (1996, and references therein). 
Based on a simple model of holes created by 
O and B stars, Oey & Clarke (1997) success- 
fully predict the observed number distribu- 
tion of SMC holes, lending support to this 
hypothesis. It should be noted, however, that 
this general picture is not without its flaws. 
Searches for the remnant stellar populations 
of the OB associations thought to be respon- 



sible for some of the holes in Ho II have not 
led to the expected result (Rhode et al. 1999). 
And in the case of the largest observed shells, 
these seem to have energy requirements sur- 
passing the output of stellar winds and su- 
pernovae. To explain those structures an al- 
ternative mechanism was proposed, the in- 
fall of gas clouds (see Tenorio-Tagle et al. 
(1987) for a numerical simulation and van der 
Hulst & Sancisi (1988) for observational evi- 
dence in the case of one of the largest holes 
in M 101). And recently, even Gamma Ray 
Bursters (GRBs) have been suggested 
possible explanation for the largest H I holes 
(Efremov, Elmegreen & Hodge 1998, Loeb & 
Perna 1998). 

The idea that supernovae deposit large 
amounts of energy into the ISM is not new. 
Cox & Smith (1974) already discussed the 
idea that the energy released by supernovae 
sets up tunnels of hot gas in the ISM. This 
picture was a substantial modification of the 
Field, Goldsmith & Habing (1969) model which 
consists of a two phase medium. The ideas 
proposed by Cox & Smith were incorporated 
within the model elaborated by McKee & Os- 
triker (1977) who suggested that the ISM has 
three distinct phases, a cool (T^ 10^ K), a 
warm (T^ 10^ K), and a hot phase (T^ 
10® K), all in pressure equilibrium. Various 
refinements have been proposed but this pic- 
ture basically still stands. 

Knowing how much kinetic energy is trans- 
ferred in expanding shells and filaments is 
of direct importance for assessing the energy 
balance in the ISM. In addition to an ana- 
lytical approach, various authors are pushing 
forward numerical simulations (e.g., Palous, 
Franco & Tenorio-Tagle 1990, Silich et al. 
1996, Palous, Tenorio-Tagle & Franco 1994), 
in order to model the evolution of the shape 
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of the H I shells in various environments, the 
fragmentation of matter on the rim of the 
shells and to investigate induced star forma- 
tion. This is important not only for tracking 
down the mechanisms which produce the H I 
shells and filaments, but also for determining 
what fraction of the ISM breaks out into the 
halo and might eventually rain back onto the 
disk. 

Our aim is to investigate the properties of 
the H I shells and bubbles in the ISM of galax- 
ies. Intimately linked to this are the ques- 
tions when fragmentation on the rim of the 
shells starts and what drives star formation. 
Unfortunately, the determination of the char- 
acteristics of H I shells in our Galaxy is not 
straightforward due to our unfavorable loca- 
tion within the disk. Looking at the near- 
est spiral galaxies (e.g., M 31, M 33, M 101, 
NGC 6946) does not improve matters all that 
much since the structures in their ISM are 
strongly modified by spiral density waves and 
differential rotation, both effects dramatically 
influencing the shape and formation of holes 
and shells. Therefore, we decided to direct 
our attention to nearby dwarf galaxies in- 
stead. Dwarfs are slow rotators, generally 
show solid body rotation, and lack density 
waves. This implies that once features like 
shells have formed, they won't be deformed 
by galactic shear and therefore tend to be long 
lived. Moreover, the overall gravitational po- 
tential of a dwarf is much smaller than in a 
normal spiral. The same amount of energy 
input of a star forming region therefore has 
a more pronounced impact on the overall ap- 
pearance of the ISM, as shown by PWBR92. 
Since dwarf galaxies also tend to have a much 
thicker disk compared to massive spirals (e.g., 
Staveley-Smith, Davies & Kinman 1992), the 
H I volume density of dwarf galaxies is lower, 



again facilitating the creation of large holes. 
And because of their smaller size, the proba- 
bility to suffer impacts from infalling clouds 
is lower compared to that for larger galaxies. 
Despite all these arguments in favor, only very 
few such H I studies of dwarf galaxies have 
been published. Beside our two neighbors, 
the SMC and LMC a detailed analysis was 
performed in one case only, that of the dwarf 
galaxy Holmberg II (PWBR92), member of 
the M 81 group. 

In this paper, we concentrate on IC 2574, a 
gas-rich dwarf galaxy in the same group. Ta- 
ble summarizes some general information on 
this object. IC 2574 (= UGC 5666 = DDO 81, 
also sometimes referred to as Coddington's 
Nebula) has been studied by several authors 
in the past. Observations with single dish 
telescopes in the 21-cm line of neutral hydro- 
gen (H I) were published by Rots (1980) and 
Huchtmeier & Richter (1988) and a first study 
using an interferometer was made by Seielstad 
& Wright (1973). A recent H I study, concen- 
trating on the dynamics and contribution to 
the total mass by dark matter was performed 
by Martimbeau, Carignan Sz Roy (1994, here- 
after abbreviated MCR94), who observed the 
galaxy with the Westerbork Synthesis Radio 
Telescope. Recent broadband optical stud- 
ies were published by Tikhonov et al. (1991). 
Studies of H II regions in IC 2574 have been 
presented by Miller & Hodge (1994), Ken- 
nicutt (1988) and Arsenault & Roy (1986). 
Spectroscopy of some prominent H II regions 
was performed by Miller & Hodge (1996) and 
Ha velocity fields were recently obtained by 
Tomita et al. (1998). In addition, four Wolf- 
Rayet stars were discovered by Drissen, Roy 
& Moffat (1993). 

We present H I observations made with the 
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NRAO0 Very Large Array (VLA) in its B- 
, C- & D-array configurations supplemented 
by optical observations obtained at the Calar 
Alto Observatory. The VLA data were origi- 
nally obtained by Puche, Brinks and Westp- 
fahl with the intention of following up on their 
work on Ho II. So far, only a low resolution 
integrated H I map was published by Puche 
& Westpfahl (1994). This paper presents the 
VLA data in a complete form and compares 
them with optical observations. The origi- 
nal data were recalibrated by us to produce 
far superior maps using the latest calibration 
and imaging techniques within AlPS. The 
data reduction and results are discussed in 
Sees. & §|. In Sec. ^, a catalog of H I holes 
and shells is presented. Sec. || discusses the 
relation between H I and Ha features. A gen- 
eral discussion is given in Sec. || and a sum- 
mary of the results is presented in Sec. |^. 

2. Observations 

2.1. Radio Observations 

IC 2574 was observed with the VLA in 
B-, C- and D-configuration. Part of the 
D-array observations were affected by solar 
interference, which was removed by delet- 
ing uu-spacings shorter than typically 0.4 kA 
and editing out bad time-ranges. In to- 
tal, 17 hours were spent on source which 
were divided into 4 hours for D-array, 2.5 
hours in C-array and 10.5 hours in B-array 
(see Table || for a detailed description of 
the VLA observations). The longer observ- 
ing time in B-array only partly compensates 
for the lower surface brightness sensitivity in 

^The National Radio Astronomy Observatory (NRAO) 
is operated by Associated Universities, Inc., un- 
der cooperative agreement with the National Science 
Foundation. 



this extended configuration. The absolute 
flux calibration was determined by observing 
1328+307 (3C286) for approximately 20 min- 
utes during each observing run, assuming a 
flux density of 14.73 Jy according to the Baars 
et al. (1977) scale. This calibrator was also 
used to derive the complex bandpass correc- 
tions. The nearby calibrators 1031+567 and 
0945+664 were used as secondary amplitude 
and phase calibrators and their fluxes were 
determined to be 1.80 Jy and 2.22 Jy, respec- 
tively. Since the systemic velocity of IC 2574 
(^ 53 kms~^) is close in velocity to H I emis- 
sion from our Galaxy, each calibrator obser- 
vation was split in two parts to form a pair 
of observations, one part having the velocity 
shifted by +300 kms~^ and the other by - 
300 kms~^ to avoid contamination by Galac- 
tic emission. In the course of the calibration 
these observations were then averaged to give 
interpolated amplitude and phase corrections. 
IC 2574 was observed using a 1.56 MHz band- 
width centered at a heliocentric velocity of 38 
kms~^. This band was divided into 128 chan- 
nels resulting in a velocity resolution of 2.58 
kms~^ after online Banning smoothing. 

The data for each array were edited and 
calibrated separately with the AlPS package^. 
The liv-data were inspected and bad data 
points due to either interference or cross talk 
between antennae were removed, after which 
the data were calibrated. We Fourier trans- 
formed our B-, C- and D-array observa- 
tions separately to assess their quality; subse- 
quently we combined all data to form a single 
dataset which was used for mapping. All the 
results presented in this paper are based on 
this combined dataset. 

Channels with velocities lower then —40 



■^The Astronomical Image Processing System (aips) 
has been developed by the NRAO. 
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kms~^ and higher than +136 kms~^ were 
found to be hne free; these channels deter- 
mined the continuum emission which was sub- 
tracted in the uv-plane. The same channels 
were used to produce a radio continuum map 
of IC 2574. 

After that, two sets of datacubes (1024 x 
1024 pixels x 68 channels each) were pro- 
duced using the task IMAGR in AlPS, each 
of them CLEANed to a level of two times 
the rms noise (Hogbom 1974, Clark 1980). 
The first data set was made with natural 
weighting, leading to a resolution of 12.4" x 
11.9" and emphasizing large scale structures. 
A second cube was produced using the RO- 
BUST weighting scheme (Briggs 1995). This 
scheme achieves the high sensitivity of natu- 
ral weighting combined with a well behaved 
synthesized beam, at a resolution close to 
that of uniform weighting. To obtain the 
optimum ROBUST value for our observation, 
the ROBUST parameter space was searched 
by producing maps in the sensitive regime of 
—2 <ROBUST< 2 and measuring for each map 
the rms noise and beamsize. Eventually, we 
chose a value of ROBUST=0.25, resulting in 
a beamsize of 6.4" x 5.9" and an rms noise 
of 0.7 mJy beam~^. A comparison with the 
beamsize of the high resolution uniform map 
(4.7" X 4.3") and the rms noise of the low- 
noise, natural map (~ 10% increase in noise) 
shows the strength of this mapping algorithm. 
Unless otherwise mentioned, the results pre- 
sented in the following are based on this RO- 
BUST cube. 

Although IC 2574 is situated at a galactic 
latitude of 43.6°, some of the channels were 
clearly confused by Galactic emission (in par- 
ticular between +10 and —10 kms~^). Since 
galactic emission shows its presence on pre- 
dominantly large scales, we were able to re- 



move the most prominent emission features 
by blanking the uv-data out to 0.1 kA in these 
channels. 

To separate real emission from the noise, 
the following procedure was applied: the nat- 
ural weighted data cube was convolved to a 
circular beam with a FWHM of 45". The 
smoothed map was then tested at the 2a level; 
if a pixel fell below this level, the counter- 
part in the cube was blanked. After that, 
the remaining peaks were inspected. Emis- 
sion that was present in 3 consecutive chan- 
nels was considered to be real while all other 
remaining spikes were considered to be noise 
and blanked. The final result was named the 
master cube. This master cube can be used as 
a mask to blank the natural and robust data 
cubes. This method ensures that the same re- 
gions are included when inspecting cubes at 
different resolutions and with different signal- 
to-noise ratios. 

2.2. Optical Observations 

IC 2574 was observed with the 2.2-m tele- 
scope of the Calar Alto Observatory^ on 1999 
January 25 and 26 in Johnson R-band (20 
min) as well as in Ha (45 min). A 2048 x 2048 
pixel CCD was used leading to 0.5" per pixel. 
The seeing was around 1.3". A focal reducer 
(CAFOS) allowed for a field of view of ra- 
dius 8' to be imaged. The usual calibration 
steps were followed using the IRAF package]^. 
Each exposure was corrected for zero offset 
(bias) and was flatfielded using skyfiats. Af- 

^The Calar Alto Observatory is operated by the Max- 
Planck-Institute for Astronomy, Heidelberg, jointly 
with the Spanish Commission for Astronomy. 

^IRAF is distributed by National Optical Astronomy 
Observatories, which is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc., 
under contract with the National Science Foundation. 
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ter that, frames were inspected for bad pix- 
els and shifted to get them ahgned. Lastly, 
the exposures were combined to eliminate cos- 
mic ray hits and to create a final image. As 
it turned out, conditions were not photomet- 
ric during the observation, hence no calibra- 
tion using standard stars was performed. In 
order to get an approximate calibration for 
IC 2574, values were taken from the litera- 
ture (MCR94). The deep R-band image was 
used for the continuum subtraction of the Ha 
map. The continuum level was determined 
by simply scaling the continuum image rel- 
ative to the line plus continuum image such 
that the foreground stars disappeared when 
the two images were subtracted. For the cal- 
ibration we relied on the Ha fluxes published 
by Miller & Hodge (1994, hereafter referred 
to as MH94). In their paper, these authors 
identify 289 H II regions in IC 2574 (for a 
thorough discussion on the Ha emission and 
its calibration see the Appendix in MH94). 
Figs. |l| and ^ show our R-band and the Ha 
image, respectively. 

Coordinates for the maps were obtained 
using the task KOORDS within the karma 
software package (for a description of karma 
see Sec. [4.11 ). We used as input an image from 
the Digital Sky Survey. Using 14 stars found 
in our CCD image, it was possible to deter- 
mine the coordinate system to an accuracy 
of better than one arcsecond, more than ade- 
quate for our needs. 

3. HI Studies 

3.1. Global H I Properties 

The channel maps obtained by the robust 
weighted cube are presented in Fig. ^ As 
can be seen, most traces of contamination by 
Galactic emission were successfully removed. 



To save space, only every other channel over 
the frequency range where line emission is 
present is shown. The heliocentric radial ve- 
locities (in kms"^) are plotted in the upper 
right-hand corner of each plane. The beam- 
size (6.4" X 5.9") is indicated in the top left- 
hand panel (note that the beamsize is so small 
that it cannot be properly reproduced in our 
maps). Many holes and shell-like structures 
are visible in the channel maps; these features 
will be the main topic of this paper. 

Obtaining the H I flux of each channel and 
hence the global H I profile of IC 2574 is 
not straightforward. As first pointed out by 
Jorsater & van Moorsel (1995), fully cleaned 
maps do not exist and any cleaned map con- 
sists of the sum of two maps: one containing 
the restored clean components and the other 
the residual map. In the former, the unit is 
Jansky per clean beam area, and in the latter 
Jansky per dirty beam area. Usually, fluxes 
are determined on a cleaned map, assuming 
that the clean beam is the correct one for the 
entire map. In reality, the flux is calculated 
correctly only for the cleaned fraction of the 
map; for the residual map the flux is overes- 
timated when the dirty beam area is bigger 
than the clean beam area (which is usually 
the case) . This becomes a real problem for ex- 
tended objects, especially in the case of com- 
bined array data, like for our VLA data set. 
For a full discussion on this topic see Jorsater 
& van Moorsel (1995). Following their pre- 
scription, the real flux of a channel map is 
given by: 

D-R 

where C is the cleaned flux, and R and D 
are the 'erroneous' residual flux and the 'er- 
roneous' dirty flux over the same area of a 
channel map, respectively. We followed this 
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procedure for every single channel to obtain 
the global H I profile as presented in Fig. ^ 
Without this correction, we would have over- 
estimated the H I flux by a factor of about 
two! 

Figure ^ shows the flux- and primary- 
beam corrected fluxes for every channel of 
IC 2574 (open squares) plotted together with 
the single dish profile published by Rots (1980) 
and the WSRT data by MCR94 (the arrows 
indicate the velocity range for which the data 
points were interpolated due to contamina- 
tion by galactic emission). From our data, 
we derive a total integrated H I flux of 373 
Jy kms~^ which should be compared with 
the value of 399 Jy kms^^ obtained by Rots 
(1980), who used the Greenbank 300 ft tele- 
scope to map IC 2574. This is excellent 
agreement, given the uncertainties of remov- 
ing Galactic emission and systematic errors 
in the absolute calibration. This also means 
that we hardly miss any H I flux due to miss- 
ing short spacings. Other values obtained by 
interferometric observations are: 440 ± 50 Jy 
kms~^ found by Seielstad & Wright (1973), 
using the Owens Valley Radio Interferome- 
ter and 286 Jy kms"^ obtained by MCR94 
from their Westerbork observations. As can 
be seen from Fig. ^, the new VLA observa- 
tions recover substantially more flux than the 
WSRT data. Adopting a distance of 3.2 Mpc, 
our flux translates to a total neutral hydro- 
gen mass of Mhi = (9.0 ± 0.1) x 10* Mq for 
IC 2574. 

3.2. Global Characteristics 

Fig. 1^ shows the H I surface density map 
of IC 2574 obtained by adding all channels 
containing H I emission of the robust cube, 
after blanking them with the master cube. 
Note that the beamsize is so small that it 



can hardly be properly reproduced. Note 
also that, although many hole-like structures 
and cavities are visible, speciflc features turn 
out to be less prominent in the total surface 
brightness map than in the channel maps. 
This is related to the fact that the galaxy 
is highly inclined, leading to overlapping of 
shell-like structures along the line-of-sight in 
the integrated map. For a more complete dis- 
cussion on that topic, see Sec. 

Fig. ^ shows the velocity fleld of IC 2574 
which was calculated in the usual manner by 
taking the flrst moment. Since the high reso- 
lution velocity fleld is dominated by small- 
scale motions within the ISM, the natural 
weighted cube was convolved to 30" (as in- 
dicated by the beamsize) to put more empha- 
sis on the global characteristics of the velocity 
fleld. The underlying greyscale is a linear rep- 
resentation of the H I surface brightness map 
at full resolution. Although the overall ve- 
locity fleld shows the typical behavior for a 
dwarf galaxy (solid body rotation in the cen- 
ter, flattening of the rotation curve towards 
the very edge of the galaxy), the general rota- 
tion is obviously disturbed by noncircular mo- 
tions (e.g., note the prominent band of lower 
velocities to the north east of the center of 
the galaxy, parallel to its major axis). These 
deviations are linked to holes in the neutral 
interstellar medium. 

Finally, Fig. |^ shows the second moment 
map, or velocity dispersion of IC 2574. Al- 
though some regions near shells reach val- 
ues as high as 15 kms~^, the overall velocity 
dispersion in the unperturbed parts is much 
lower. Using H I proflles, we derive the veloc- 
ity dispersion in the disk by looking at quies- 
cent regions. Several lines of sight were aver- 
aged to give a one sigma velocity dispersion 
of 7 ± lkms~^ (similar to what is found in 
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Ho II). This shows that our velocity resolution 
of 2.5 kms~^ (FWHM) was quite sufficient to 
resolve the lines of neutral hydrogen. 

In this paper, we focus on the small-scale 
structure of the ISM in IC 2574. For our 
analysis, however, we need to adopt values 
for the orientation of the object which are 
usually derived on the basis of a kinematical 
analysis (rotation curve fitting). A full anal- 
ysis of the dynamics of IC 2574 based on our 
observations, however, is beyond the scope 
of this paper. Luckily, an adequate analysis 
has already been published (see the paper by 
MCR94). In the course of their study they 
derive a systemic velocity of 58 ± 7 kms~^, 
adopt an inclination of 75° ± 7° and a posi- 
tion angle of 52°. Depending on their dynam- 
ical model, they find that 70 — 90% of the 
total mass of the galaxy is in the form of dark 
matter. Following the same procedure as per- 
formed for the dwarf galaxy II Zw 33 (Walter 
et al. 1997) we carried through a first analysis 
of our data confirming the results by MCR94 
regarding the rotation curve and the global 
orientation of IC 2574. For ease of comparison 
we therefore adopted their values throughout 
this paper. 

3.3. The scale height of the H I disk 

The scale height h of the H I layer of 
IC 2574 plays an important role in the de- 
termination of the properties of the holes, 
e.g., when translating the observed H I sur- 
face density into an H I volume density. Here 
h is defined as the la value for a Gaussian 
distribution of the H I layer. Recent stud- 
ies have shown that in general dwarf galaxies 
have thicker H I disks than more massive spi- 
ral galaxies. Their aspect ratio tends to be 
of the order of 10:1 rather than 100:1, e.g.. 
Ho II: h = 625 pc (PWBR92) or NGC 5023: 



h = 460 pc (Bottema et al. 1986; see also 
Staveley-Smith et al. 1992). In the case of 
IC 2574 there are also indications that the 
H I is distributed in a thicker disk. MCR94 
already noticed that, looking at the appar- 
ent axial ratio of the integrated H I map and 
assuming a real inclination of 75°, the H I 
must be distributed in a somewhat thicker 
disk. Moreover, the fact that we observe, 
as we will argue below, spherically expanding 
shells with diameters of up to almost 1 kpc 
implies that the extent of the H I layer must 
be of the same order of magnitude. Hence, 
the scale height of the H I gas should be of 
order 350-400 pc. 

An approximate value for the scale height 
can be obtained fairly simply, as shown by 
PWBR92. We follow their approach here. 
The scale height is proportional to the veloc- 
ity dispersion of the gas and inversely propor- 
tional to the square root of the mass density 
in the disk (Kellmann 1972; van der Kruit 
1981). For a density distribution that can be 
described as p{z,R) = p{0, R) sech'^{z/ zq), Zq 
is given by 

Zo{R) = I ^^"^ 

^27rG'p(0,i?) 

where R is the galactocentric distance, agas 
the velocity dispersion of the gas in quies- 
cent regions of the galaxy (Sec. ^.2[ ), and p{R) 
the total volume density (gas, stars and dark 
matter) of the galaxy in the plane. However, 
since seQb?{z/zQ) ^ ^^^{—z^ / z^) (to within 
5% in the region of interest) we can also as- 
sume a Gaussian distribution in z. The Icr 
scale height h is then defined as /i = zq/\/2. 
To calculate the average density of IC 2574 
in the plane we assume that the total mass 
of the galaxy can be estimated from the last 
measured point {Rmax) of the observed ro- 
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tation curve (MCR94) and that this is dis- 
tributed over a spherical volume of radius 
Rmax- MCR94 derive a rotation curve with 
v{Rmax) = 66 kms"^ at Rmax = 8.0 kpc. 
This corresponds to a total dynamical mass 
of 8.3 X 10^ Mq. Distributing this mass over 
a spherical volume with radius Rmax gives an 
average value for p of 0.0039 Mq pc~^. Since 
this is only the mean density of the galaxy 
and the density near the plane is certainly 
higher, we take twice this derived value (see 
PWBR92 for a justification) to get a first ap- 
proximation to the real volume density p{R) 
near the center. Using this approach, we find 
p{R) = 0.008 Mq pc-3. 

Substituting the observed velocity disper- 
sion of 7 kms~^ (see Sec. [HT^ ) in the equa- 
tion given above we derive a scale height of 
h = 345 pc for IC 2574. Note that this 
value, given all the uncertainties in the deriva- 
tion, is in good agreement with the value de- 
rived based on the sizes of the largest holes 
{h = 350 — 400 pc). We therefore adopted a 
scale height of h = 350 pc for the H I layer of 
IC 2574. 

We can now calculate the volume density 
of the H I from the projected and corrected to 
face-on H I surface density. The surface den- 
sity profile was obtained by elhptically aver- 
aging the flux- and primary-beam corrected 
total H I map (assuming an inclination of 
75°). Again, this is correct to first order only 
as in the case of a thick H I distribution the 
deprojection of the observed projected H I 
distribution is not trivial. We then calculated 
the volume density nHi{R) of the H I in the 
plane of the disk using 

/+00 — ^2 

^ nniiR) exp (^) dz = V27r/i 

4. H I Holes in IC 2574 



4.1. The H I— hole catalogue 

We conducted a search for H I structures 
in IC 2574 using the tasks KVIEW, KPVS- 
LICE and KSHELL that are built into the re- 
cently developed karma^ visualization soft- 
ware package. KVIEW offers many different 
ways to look at a movie of a datacube whereas 
KPVSLICE produces real-time position veloc- 
ity (pV) cuts through a cube in any orienta- 
tion which tremendously facilitates the search 
for and identification of hole-like structures. 
Finally, kshell was used to make radius- 
velocity plots of some of the more prominent 
shells and holes. 

First, the ROBUST cube was searched for 
the smaller holes in IC 2574. In order to 
be sensitive to more extended structures and 
therefore bigger holes, we smoothed our ro- 
bust cube to 11", 15" and 30", respectively, 
and inspected each of them. For each H I 
hole in IC 2574 we determined its position, 
size and expansion velocity. Both authors in- 
dependetly produced a list of holes to reduce 
personal bias, based on some simple crite- 
ria regarding the determination of the center, 
size, expansion velocity and systemic velocity 
of the holes. After that, the two lists were 
compared and merged. There was an overlap 
of about 75% in terms of well defined holes 
which lends some credibility to the identifi- 
cation process. We reexamined the remain- 
ing structures and, erring on the side of cau- 
tion, ended up with a catalogue of 48 holes. 
Although this list cannot be free from per- 
sonal bias, the method closely follows that 
employed in the case of M 31 and Ho II which 
should make a comparison between those ob- 

KARMA visualization software package was de- 
veloped by Richard Gooch of the Australia Telescope 
National Facility (ATNF) 
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jects and IC 2574 valid. 

It is interesting, at this point, to reflect 
on the reahty of the features. Do they really 
constitute large-scale, coherent structures or 
are we tricked by nature into seeing coher- 
ence in what is basically a turbulent, per- 
haps fractal medium? Mac Low et al. (1998) 
recently produced sets of model data cubes, 
similar to those obtained via radio (H I) or 
optical (Ha) observations (i.e., position, po- 
sition and velocity). We analysed, using the 
same techniques as explained above, a cube 
of uniform, isotropic, isothermal, supersonic, 
super-Alfvenic, decaying turbulence. The 
simulations assume that the gas is optically 
thin and emits with a strength directly pro- 
portional to its density. The model corre- 
sponds to a weakly magnetized ISM. Single 
cuts through the data cube show features re- 
markably similar to the ones seen in the H I 
data. However, when focussing on a potential 
feature and investigating its morphology in 
position-position or position-velocity space, 
no coherent structure emerges. This lends 
us to believe that shells due to random (or 
at least turbulent) processes make up a small 
fraction of the holes we classifled. Instead, it 
is more likely that the H I holes and shells 
are due to mechanisms which act on large 
scales, such as stellar winds and (multiple) 
super novae. 

Following Brinks & Bajaja (1986), we dis- 
tinguished between three different types of 
holes based on their appearance in the pV- 
diagrams. A Type 1 hole corresponds to a 
total blowout (i.e., we don't observe the re- 
ceding nor the approaching side of the shell). 
In that case, it is impossible to determine the 
expansion velocity and hence the kinematics 
of the hole. Figure ||c shows a pV-diagram 
of a typical Type 1 hole (our hole number 30) 



using KPVSLICE. 

A Type 2 hole corresponds to a hole which 
is offset with respect to the plane of the 
galaxy. In that case, we witness a deformation 
in the pV-diagram towards the plane of the 
galaxy and a blowout to the opposite side (see 
Silich et al. 1996 for numerical simulations of 
such a case — their case C). The amount of 
deformation yields the expansion velocity of 
the shell towards the midplane. The part of 
the shell that points away from the plane is 
usually not visible due to lack of sensitivity. 
Figure |^d shows a pV-diagram of a typical 
Type 2 hole (our hole number 33). 

Type 3 is the 'classic hole' and corresponds 
to a more or less complete expanding shell. 
The signature in pV-space is that of an ellip- 
tical hole. The expansion velocity is equal to 
half the difference between the velocity of the 
approaching and receding sides. The systemic 
velocity of the hole is deflned by its center. 
Figure ^ shows a pV-diagram of a typical ex- 
ample of a Type 3 hole (our hole number 21). 
Figure ||b shows the same hole, but analyzed 
with the task KSHELL. In this flgure, the hor- 
izontal axis corresponds to radius r (distance 
from the center of the hole) and the inten- 
sity corresponds to the average H I density 
within an annulus with radius r around the 
center of the hole. The vertical axis is again 
velocity. Since we are looking at a radius- 
velocity diagram, we will call such a repre- 
sentation an rV-diagram in the following. For 
an ideal shell, a half ellipse should be visible 
in which the diameter in the f-direction is 
twice the expansion velocity and its dimen- 
sion in the r-direction the radius of the shell. 
The advantage of looking at a near circular 
H I shell with kshell lies in a substantial in- 
crease in signal-to-noise ratio as one averages 
over many pixels whereas in a pV-diagram 
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one merely cuts through an H I structure. 

Note that in the case of a Type 1 hole, the 
absence of detectable emission corresponding 
to the front and back of a shell causes a hole 
in rV-space. For obvious reasons, a Type 2 
hole does not show a prominent signature in 
an rV-diagram either. 

As was already mentioned earlier, the holes 
were found to be much more prominent in the 
channel maps than in the integrated H I map. 
This is due to the high inclination of IC 2574 
which leads to several H I shells at different lo- 
cations within the disk and possibly different 
velocities overlapping along the line-of-sight, 
as is the case in M 31. The total H I surface 
brightness map was therefore of little use for 
the identification of a hole. Note that in the 
case of Ho II, only the integrated map was 
used because of the much lower inclination of 
the galaxy. This implies that the approach to 
search for holes needs to be adapted to the 
orientation of the galaxy. 

Figure ^ shows a plot of the distribution of 
the 48 H I holes detected with confidence in 
IC 2574. The greyscale map is a linear repre- 
sentation of the total H I surface brightness. 
Some of the holes are clearly visible in the to- 
tal H I map but the majority of holes doesn't 
show up at all and can only be detected in 
the channel maps or pV-diagrams. 

The final list of the detected H I holes is 
given in Table ^. The properties listed are as 
follows: 

Column 1: Number of the hole (increasing 
with increasing right ascension) 

Column 2, 3: Position of the center of the 
hole in right ascension and declination, 
respectively (B1950.0 coordinates). Typ- 
ical errors are 3". 



Column 4: Heliocentric velocity fhei of the 
hole defined by the channel in which the 
hole is largest or most prominent. The 
accuracy of fhei is about 2.5 kms~^ (the 
velocity resolution of our observations). 

Column 5: Diameter d of the hole in kpc. 
Note that most of the holes are circu- 
lar. In the case of elliptical holes, the 
diameter along the major axis, along 
the minor axis and the position angle 

(2 X fejnaj, 2 X brain, P^) IS givCU (but SCC 

Sec. [4.2| ). The accuracy of d is around 
95 pc (the FWHM of our synthesized 
beam) . 



Column 6: The expansion velocity fcxp 



of 



the hole in kms~^. Note that in the 
case of a Type 1 hole, no expansion ve- 
locity can be measured. Typical errors 
are about 2 kms~^. 

Column 7: The estimated H I volume den- 
sity rim (in pc~'^) of the region where 
the hole is situated before the hole was 
created, nni was calculated based on 
the H I surface density, as described 
in Sec. Given the uncertainties in 
the deprojection and the assumptions 
which went into the derivation of the 
mass density within the disk, we esti- 
mate a relative error in the value of nm 
to be of order 30%. 

Column 8: The type of the hole (1-3), as 



defined in Sec. O. 



Table ^ summarizes the physical properties 
derived from Table |^. The properties listed 
are as follows: 

Column 1: Number of the hole 
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Column 2: Diameter d of the hole. For el- 
liptical holes, the diameter was defined 
using the geometric mean 



d 



2 &maj X 2 6j] 



Column 3: Expansion velocity fexp of the 
hole in kms~^. Same as column 6 in 
Table |. 

Column 4: Kinematic age of the hole in units 
of 10^ years. This age determination as- 
sumes that over its entire lifetime the 
hole was expanding with v^xp- Hence: 



t 



id/2) 



^exp 



The accuracy of the kinetic age is typi- 
cally 20%. Because of the assumption of 
a constant expansion velocity, the kine- 
matic age is an upper limit. 

Column 5: Indicative H I mass Mhi in units 
of 10^ M0. It is the mass which filled 
the observed hole under the following 
assumptions: i) the hole is spherical, ii) 
the gas was distributed uniformly in the 
galaxy, and iii) the hole is completely 
empty now. Note that if one can assume 
that only a minor fraction of the sur- 
rounding H I is ionized, this mass cor- 
responds also to the total H I mass on 
the rim of the hole. Mhi was calculated 
using 



M, 



ni = -7T{d/2fnm{R) 



where n{R) is the volume density of 
the gas at the galactocentric distance R 
(column 7 in Table To correct this 
value for the amount of primordial he- 
lium, a factor of 1.33 has to be applied. 



The relative error is, as d enters to the 
third power, of order 30%. Note that in 
the case of the largest holes, this value 
is more likely to be an upper limit since 
TT-Hi decreases with increasing distance z 
from the plane. 

Column 6: The initial total energy Ee that 
was released by the supernova explo- 
sions which created the hole (in units 
of 10^*^ ergs). Note that the kinetic en- 
ergy of the expanding shell is about 1% 
to 10% of Ee (depending on the am- 
bient density, see Chevalier 1974). To 
determine E^, we used the equation by 
Chevalier (1974) who derived a fit to 
numerical results of an expanding shell 
created by a single explosion (see also 
McCray & Snow 1979). Expressing E^ 
in terms of observable quantities d, fexp 
and uq, he obtained 



E. 



5.3 X 10^3^1-12 {d/2f'^v, 



3.12 „, 1.4 

exp 



ergs 



where fexp is the expansion velocity in 
kms~^, no the density of the ambient 
medium in particles per cubic centime- 
ter, and d is the diameter of the shell in 
parsecs. In our calculations we replaced 
no by nni as was done, e.g., for M 31 
and Ho II. This is not strictly correct 
as we are ignoring contributions by He 
and heavier elements. To correct E-^ for 
this contribution, all values should be 
multiplied by a factor of ~ 1.5. In addi- 
tion, the derived H I densities are aver- 
age values for a given galactocentric dis- 
tance which introduces another source 
of uncertainty. Note that the numeri- 
cal coefficient in Chevaliers equation de- 
pends on the radiative energy loss rate 
and thus on the heavy element abun- 
dance in the expanding shell. Not only 
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is this abundance different for various 
galaxies, but it decreases with galac- 
tic radius, both for our own and exter- 
nal galaxies and by more than an order 
of magnitude (e.g., Heiles 1979). The 
last phase in the expansion of a shell is 
reached when, after the shell has con- 
tinued to collect more of the ambient 
medium it will slow down to velocities 
comparable to the random motions of 
the ISM (around 7 kms~^) and hence 
stops expanding. Most probably the 
largest structures we observe are in this 
latter phase (see Sec. [4.3[ ). From this 



brief discussion it should be clear that 
the derived energies Ee are only order 
of magnitude estimates. 

4.2. Notes on single holes 

The following holes deserve a more detailed 
description: 

Hole 7^35: This hole causes a major distur- 
bance in the pV-diagrams of IC 2574. 
It is one of the few elliptical holes, with 
dimensions of 1000 x 500 pc. Although 
there is some circumstantial evidence 
that it is a superposition of two spheri- 
cal holes, it was treated as one hole. As 
we will show in Sec. |, the most promi- 
nent H II region in IC 2574 is located 
on the rim of this supergiant shell. The 
shell is filled with soft X-ray emission as 
observed with the ROSAT telescope. A 
detailed multiwavelength analysis (opti- 
cal, radio. X-ray) of this particular su- 
pergiant shell has been presented else- 
where (Walter et al. 1998) 

Hole #22: Catalogued elliptical in shape, this 
hole is actually a superposition of sev- 
eral holes, leading to the overall impres- 



sion of an elliptical feature; we treated 
this as a single object since subdividing 
it in, say, three holes would have been 
arbitrary. 

Holes 7^ 13, 23, 26: As can be seen from 
the channel maps (Fig. ^), some holes 
are not clear-cut cases in terms of the 



definition of a hole given in Sec. 4.1 



They are possibly chance superposition 
of old shells and fragments leading to 
the impression of a hole in the inte- 
grated H I map (see also Sec. ^^ ). We 
therefore decided not to include these 
holes in the statistical analysis. 

4.3. Notes on large— scale structures 

As can be seen from Table 0, only few 
shells larger then 1000 pc have been cat- 
alogued in IC 2574. It has to be mentioned, 
however, that there appear coherent features 
larger then 1 kpc in the channel maps. Several 
of these can be identified in the channel maps 
(Fig. ^), e.g., between 45 and 70 kms~^. Note 
that these structures are related to the major 
disturbance in the velocity field due north of 
the center of the galaxy. Obviously, it is very 
difficult to catalogue these features as they 
lack the clear coherence, both in position as 
well as in velocity, of H I shells. Moreover, 
because of their large size, their dynamics is 
not dominated by some leftover expansion but 
rather by large scale galactic rotation modi- 
fied by the effects of nearby, smaller holes. 
The reason why we wish to mention these 
large scale features is because along some of 
the apparent arcs we find H II regions. 

5. The H I-Ha connection in IC 2574 
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5.1. Global correlations 

Until now, we have focussed on the struc- 
ture and dynamics of the neutral ISM of 
IC 2574. In the following, we will compare 
these results with our Ha observations which 
reveal the sites of recent or ongoing SF in 
IC 2574. As noted earlier (Sec. P, the ex- 
panding H I structures are widely thought to 
result from recent star formation. They sweep 
up gas causing the volume density on the rim 
of the shell to increase until gravitational in- 
stabilities lead to molecular cloud formation 
and hence secondary SF (e.g., Palous et al. 
1990, Elmegreen 1994). Therefore, a compar- 
ison of the position of H I holes with H II 
regions should show some correlation, the Ha 
emission being predominantly found along the 
rims of larger holes, as evidence for propagat- 
ing SF. In the case of the smaller, young holes 
one would expect some to be still filled with 
diffuse Ha emission coming from the ionized 
wind-blown cavities around the still present, 
young O and B stars. 



Figure [T^ shows an overlay of the Ha map 
with the H I holes we detected in our sur- 
vey (Sec. |4.1j ). The width of the rims of the 
H I holes is of course in reality larger then 
indicated in this plot. The subregions A-C 
will be discussed in more detail in the follow- 
ing section. As can be seen from this overlay, 
there is a prominent correlation of ongoing SF 
with the rims of the large H I shells which is 
in general agreement with the picture painted 
above. There is a strong tendency for H II re- 
gions to be situated on the rims of the H I 
shells (for some very nice examples see holes 
8 and 35) and only in a very few cases diffuse 
Ha emission does not seem to be related to 
an H I structure. Note that in one case only 
(hole 32), diffuse Ha emission (partly) fills an 
H I hole. 



Fig. 11 presents an overlay of the Ha map 
with one contour of the H I surface bright- 
ness map (flux and primary beam corrected) 
plotted at 1.7 x 10^^ atoms cm~^. As can be 
seen from this overlay. Ha emission is globally 
correlated with a high H I surface brightness 
(for an exception in the very north see below). 
Kennicutt (1989) used the disk instability re- 
lation of Toomre (1964) to show that, for spi- 
rals, there is a threshold density for star for- 
mation. For gas densities above the threshold, 
the star formation rate (SFR) is consistent 
with a Schmidt law. Near the threshold, the 
correlation breaks down, and one sees bursts 
of SF. Much below the threshold the SFR is 
extremely low. Skillman (1987) has observa- 
tionally determined this threshold for dwarf 
galaxies to be on the order of 1 x 10^^ atoms 
cm~^ (at a resolution of 500 pc) which is in 
general agreement with our results (see also 
below). For a detailed discussion on the mod- 
els that describe the relationship between gas, 
stars and star formation in irregular galaxies, 
the reader is referred to the paper by Hunter, 
Elmegreen & Baker (1998). 

5.2. Analysis of individual H II regions 

Using our high resolution H I data of 
IC 2574 we are now able to perform a much 
more detailed analysis and comparison of the 
properties of the H II regions with the H I fea- 
tures that are related to them. When making 
overlays of the Ha map with the high resolu- 
tion cube we found that in virtually all cases, 
single H II regions are correlated with distinct 
H I structures like clumps or clouds. We de- 
fined 40 regions where Ha emission and H I 
clouds seem to be related. These regions are 



outlined by the squares in Figs. p!4-[14| which 
are meant to indicate which H II regions were 
lumped together. They are numbered 1' to 
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40' to avoid confusion with the numbers of the 
H I holes. Note, that sometimes two or more 
Ha features belong to the same H I struc- 
ture (like a clump) and were thus considered 
as one region (e.g., in region 40'), whereas 
in other cases individual H II regions were 
treated as separate entities (e.g., regions 11' 
and 12'). The Ha fluxes of the regions tell 
us something about the current SF rate and 
the average number density of the ionized gas 
whereas the H I data yield approximations 
for the number density of neutral hydrogen, 
the velocity dispersion of the gas in that re- 



gion and the thresholds for SF. Figs. p!^-P 
are enlargements of the regions A-C marked 
in Fig. ^ and display the regions defined by 
us. The position and size of the H I holes are 
plotted as well for comparison (please refer to 
Fig. for cross referencing of the holes). The 
fluxes of the H II regions were taken from the 
observations of MH94. Table ^ summarizes 
the cross identifications between their and our 
regions. In this table we list the coordinates 
of the centers of our regions as well. 

Table ^ summarizes the results of our anal- 
ysis of the H II regions. The first column 
refers to the number of the H II region as 
defined in Figs. |T2|-|Ill From the fiuxes pub- 
lished by MH94, we deduced the total fiux 
F{Ha) of our groups by simply summing 
their fiuxes in the respective regions (col- 
umn 2, Table |]). From this number, the Ha 
luminosity can be calculated using L{Ha) = 
4:TiD^F{Ha), which, in the case of IC 2574 
(^=3.2 Mpc), translates to L{Ha) = 3.1 x 
lO^^F(Ha) where F{Ha) is given in ergs cm~^ s 
and L{Ha) in units of solar luminosities Lq 
{ILq = 3.85 X 10^^ ergs s~^). The num- 
ber of Lyman continuum photons per second 
is given by iVLyc = 2.82 x 10^^ L(Ea). Ta- 
ble ^ (see Devereux et al. 1997 and references 



therein) summarizes the expected properties 
of high mass stars; column 1 gives the spec- 
tral type, column 2 the stellar mass, column 3 
the bolometric stellar luminosity, column 4 
the number of Lyman continuum photons per 
second, column 5 the Ha luminosity L(ifa), 
and column 6 the ratio of bolometric to Ha 
luminosity. These equations and values can 
be used to estimate, for example, the equiva- 
lent number of 05 stars, for each region (see 
below). 

Column 3 in Table ^ presents the volume 
density of the ionized free electrons over the 
area where we find Ha emission (as defined 
by MH94) . This was arrived at by first calcu- 
lating the emission measure (EM) of each of 
their regions according to: 

EM = 2.41 X 10^ r°-^2 F{Ha) cm"^ pc 

(Peimbert et al. 1975), where F{Ha) is the 
Ha fiux in units of ergs cm~^ s^^ sterad^^ 
and T is assumed to be the canonical value of 
10"^ K. From that, the mean electron density 
for every single region published by MH94 was 
derived using 



EM 



pc 



where Ue is the electron density and d the 
path length in pc through the H II region 
(assuming a spherical distribution of the Ha 
emission). Since the derived values for 
were very similar for a single region, we cal- 
culated the value for the electron density for 
each of our regions by taking an intensity- 
weighted average of the contributing MH94 
regions. The average electron densities range 
from 0.5 to 2 cm~^. This is, most likely, a 
lower limit for the actual electron density due 
to internal absorption within IC 2574. Note 
that in the above, effects due to the filling fac- 
tor of the gas are ignored. The real densities 
in the H II regions are likely much higher. 
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In order to derive properties of the H I fea- 
tures that are physically related to the H II 
regions, we defined areas A around each of 
the H II regions with an effective radius given 

by reff - 



A /it, where is given in pc (col- 
umn 4). Although this is a somewhat subjec- 
tive process, we found that slight variations 
in the size of the aperture didn't change the 
derived properties by more than 20 — 30%. 
Column 5 represents the total H I mass found 
within the aperture using the flux and pri- 
mary beam corrected H I data. For every 
single region, we measured the velocity dis- 
persion of the H I by taking spectra through 
the H I cube. Column 6 summarizes the mea- 
sured velocity widths defined as the full width 
at half maximum (FWHM) from which the 
real velocity dispersion (assuming a Gaus- 
sian velocity distribution) can be calculated 
using Av = 2 a/2 In 2 cTv = 2.355 dv Interest- 
ingly, these values tend to be high where two 
or more shells overlap (e.g., the H II com- 
plexes 16, 31, 32, 33). These high disper- 
sions are probably caused by blending along 
the line of sight of these shells. 

Column 7 presents the average density of 
the neutral hydrogen which was calculated us- 
ing ^ 

PHI = MHi/(-7rr3g) 

with Mhi and les the values from columns 5 
and 4, respectively. Values for pni range from 
about 1 cm"^ to values as high as 6 cm~^. 
Note that these are only averages and that 
locally the H I density might be more en- 
hanced. Based on the analysis by MCR94 and 
some crude mass modeling (Sec. |0| ) we de- 
rive an average density in the disk of 0.008 
MqPc"'^ = 0.3 atom cm~^. The regions in 
which star formation is currently active thus 
have an overdensity of an order of magnitude 



(see Table |]). 

The average surface density of the H I re- 
gions is shown in column 8. Note that virtu- 
ally all values are larger than 1 x 10^^ cm~^ 
which is in excellent agreement with the ob- 
servational value given by Skillman (1987). In 
other words, SF is taking place if the pro- 
jected H I surface density reaches values of 
about 10^^ cm~^ not only for a resolution of 
500 pc, as observed by Skillman but also for 
almost an order of magnitude improvement 
(~ 95 pc)! The highest surface density found 
is of order 2.5 x 10^^ cm~^. 

Finally Column 9 summarizes some re- 
marks regarding the regions. Most of the 
H II regions are related to the rim of an H I 
shell. It should be noted at this point that 
in the course of checking for correlations be- 
tween Ha and H I emission, we detected 5 
new tiny holes at the positions were there is 
Ha emission. However, as these holes were 
not detected in our inspection of the H I cubes 
we didn't include them in the final catalogue 
which is thought to be complete for holes with 
radii larger then about 100 pc only (as dis- 



cussed in Sees. 4.1 andBT 



We catalogued four of the H I features as 
"clumps" (regions 1', 2', 6' & 22'). By 'clump' 
we mean a well-defined H I region in 2 spa- 
tial dimensions as well as in velocity space. 
These clumps look like small, filled circles 
in position-velocity space. For those cases 
we asked ourselves whether these clumps are 
gravitationally confined or not. To answer 
that question, we calculated the virial masses 
of these objects using 

M^ir = 250 Av^ reff Mq 

(Rohlfs & Wilson 1996), where Av is the full 
velocity width at half maximum (FWHM) in 
km s~^ and Teg the radius of the clump in 
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pc. The virial masses of the four regions 
are: M^i,{V) = 13 x IO^Mq, M^i,{2') = 
38 X IO^Mq, Mvir(6') = 21 X IO^Mq and 
M™(22') = 7 X IO^Mq. As we see, the virial 
masses are all about one order of magnitude 
higher than the H I masses. So, the clumps 
are either not gravitationally bound or they 
are held together by as yet undetected ma- 
terial. One possibility is that they contain 
an order of magnitude more mass in the form 
of molecules. This will be hard to confirm, 
however, as dwarf galaxies have been notori- 
ously difficult to detect in CO, the most abun- 
dant tracer for molecular gas (e.g.. Sage et 
al. 1992). Dark Matter could play a role, 
but the mass density, including DM, near 
the above mentioned SF regions is of order 
0.008 MqPc~^ (see Sec. ^.3[ ), insufficient by a 
wide margin. 

If the observed velocity dispersion is not 
due to gravitation, the question remains, why 
do we see clumps? We can think of two more 
mechanisms, magnetic confinement or pres- 
sure. The former is unlikely to play a large 
role. Dwarf galaxies do not possess large 
scale magnetic fields (the exceptions being 
NGC 4449 and IC 10, e.g., Chyzy et al. 1998), 
any non-thermal radio continuum radiation 
being confined to the strongest regions of ac- 
tive star formation (Deeg et al. 1993). 

We did not find any correlation of Ha fiux 
with the surrounding H I density. However, 
one wouldn't really expect any strong corre- 
lation of this kind. Not only is the Ha emis- 
sivity of a star forming region strongly depen- 
dent on the age of the region (see, e.g., Ger- 
ritsen & Icke 1997), but it is also dangerous 
to try to deduce actual physical properties of 
the SF regions from mean values, like average 
electron or H I densities, while ignoring spa- 
tial variations in their volume filling factors. 



In some regions within IC 2574 (e.g., re- 
gion 40'), we find high H I surface densi- 
ties but only little Ha emission indicating 
that possibly a major star formation burst 
is about to commence there and that there- 
fore most of the gas is still neutral. With 
the help of the equations given above we 
can calculate how many stars are approxi- 
mately needed to create the observed Ha fiux. 
In the case of region 40' we find a fiux of 
1.5 X 10"^*^ ergs cm~^ s^^ and deduce a lumi- 
nosity of 4.7 X 10^ Lq which is the equivalent 
of a few 07 stars only (see Table Note, 
however, that there exist also regions where 
the situation is reversed and most of the H I is 
being ionized by the Lyman continuum pho- 
tons of the most massive stars of the power- 
ful H II regions. Region 34' is a prominent 



example (see also Fig. |TTD. From the Ha fiux 



of F{Ha) 1 X 10~ ergs cm~ s~ we derive 
an Ha luminosity of L{Ha) = 3.1 x 10^ 
which corresponds to an equivalent of 22 05 
stars, in other words a major SF region. 

6. Discussion 

6.1. What created the H I holes? 

The ISM of IC 2574 is dominated by ex- 
panding H I holes and shells. Similar struc- 
tures have been found in other gas-rich dwarf 
galaxies (e.g., Holmberg II (PWBR92), IC 10 
(Wilcots & Miller 1998) and in more mas- 
sive spiral galaxies like our Galaxy (Heiles 
1979), M 31 (Brinks & Bajaja 1988) and M 33 
(Deul & den Hartog 1990). One of the major 
questions concerning these structures is what 
physical process or processes are responsible 
for their origin. 

Heiles (1979), who was the first to rec- 
ognize expanding H I shells in the Galaxy, 
concluded that 'many of the H I structures 
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[...] were probably produced by large amounts 

of energy injected directly to the interstellar 
medium'. One obvious mechanism is the after 
effects of massive star formation: strong stel- 
lar winds and supernova explosions (see the 
review by Tenorio-Tagle and Bodenheimer 
1988). Detailed numerical modeling by var- 
ious authors (e.g., Palous, Franco & Tenorio- 
Tagle 1990, Silich et al. 1996, Palous, Tenorio- 
Tagle & Franco 1994) has strengthened con- 
fidence in this explanation, as has the work 
by Oey & Clarke (1997) who in a (semi- 
) analytical way managed to predict the H I 
shell size distribution in the SMC based on 
the observed Ha luminosity function. 

Having said that, there are potentially prob- 
lems with this interpretation, as mentioned 
by Radice, Salzer & Westpfahl (1995) and 
Rhode, Salzer & Westpfahl (1999). Looking 
for remnant stellar clusters at the centers of 
several large holes in Holmberg H they failed 
to find a population of A and F main sequence 
stars. Moreover, it is difficult to imagine that 
the largest H I shells which are located be- 
yond the optical galaxy should have been cre- 
ated by star formation (an argument which 
holds equally well for some of the holes in 
IC 2574). This is corroborated by ultravio- 
let imaging of Holmberg 11 by Stewart (1999) 
who finds no FUV emission from the outer 
regions of Holmberg H. 

It should be noted, though, that the non- 
detection of population A and F stars does 
not necessarily mean that none of the H I 
are created by SF. In dwarf galaxies, in the 
absence of shear, H I shells can reach large 
dimensions and grow so old that no obvious 
bright cluster remnant will remain visible in 
the center (Efremov, Elmegreen & Hodge ). 
Some of the observed supershells are so old 
that the brightest members of stellar cluster 



which created the hole may have evolved off 
the main sequence and dispersed (assuming 
a similar IMF as in the Galaxy). This may 
make a detection of old clusters difficult. Such 
a case (LMC 4/Constellation HI) is described 
by Efremov & Elmegreen (1998). 

Another potential problem with the con- 
ventional picture for the origin of the holes is 
that for the largest ones the energy require- 
ments seem to surpass the amount that can 
possibly be released into the ambient ISM via 
stellar winds and supernova explosions. The 
giant shell in M 101 (van der Hulst & Sancisi 
1988) is just one example of such a structure. 
Tenorio-Tagle (1987) proposed as an alterna- 
tive the infall of high velocity clouds which 
could deposit an amount of energy of 10^^ 
to 10^^ ergs per collision (see the review by 
Tenorio-Tagle & Bodenheimer 1988 for fur- 
ther references) . Even though this mechanism 
seems a viable explanation for at least some 
holes, it is unlikely that it plays a dominant 
role. Certainly it can't be as important in 
dwarf galaxies because of their much smaller 
collisional cross section. 

Recently, an alternative explanation for 
the creation of supergiant H I shells in galax- 
ies has been proposed, namely GRBs or Gamma 
Ray Bursters (Efremov, Elmegreen & Hodge 
1998, Loeb & Perna 1998). Even though the 
basic physics of GRB are not yet known the 
picture to emerge is that they are highly en- 
ergetic events (of order 10^^ ergs) that release 
huge amounts of energy in a very short period 
of time (order of a few seconds). 

Although we recognize that perhaps not all 
holes can be explained by stellar winds and 
supernova explosions, we feel that there is suf- 
ficient evidence for the conventional interpre- 
tation to warrant an analysis in those terms, 
which is what we pretend in the rest of this 
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paper. 

6.2. Statistical properties of the holes: 
a comparison with other galaxies 

For the first time we are now in a position 
to address an important question regarding 
the creation and evolution of hole-like struc- 
tures in the interstellar medium: how are 
their properties related to their environment, 
i.e., the type of host galaxy. In what follows 
we will compare the observed and derived H I 
hole properties of IC 2574 with those found 
in M 31 (Brinks & Bajaja 1986), an exam- 
ple of a massive spiral galaxy similar to our 
own, M 33 (Deul & den Hartog 1990), a less 
massive spiral and Ho 11 (PWBR92), another 
dwarf galaxy in the same group of galaxies, 
yet four times less massive than IC 2574. In 
other words, the sequence M 31 - M 33 - 
IC 2574 - Ho II spans a large range of dif- 
ferent Hubble types from massive spirals to 
low-mass dwarfs. 

Despite the fact that a similar study, and 
at considerably higher sensitivity, is now avail- 
able for the SMC (Staveley-Smith et al. 1997), 
we decided not to include their results. The 
main reason being that the linear scales (or 
spatial frequencies) sampled by the ATNF 
hardly overlap with those observed in the 
galaxies listed above. The linear resolution, 
at 28 pc, is almost four times higher than, 
for example, the VLA maps on IC 2574. At 
the other end of the spectrum, because of the 
lack of short spacing information, structures 
larger than a few hundred parsec will have 
been missed. Moreover, the SMC is a very 
disturbed system, being torn apart by tidal 
forces due to interactions with the LMC and 
the Galaxy. 

Limiting the comparison to the four galax- 
ies listed above has some further advantages. 



The linear resolutions are very similar, as are 
the velocity resolutions with which they have 
been observed (see Tab. ?? for a summary). 
In addition, all four galaxies were examined 
in more or less the same fashion, one of the 
authors (EB) having taken part in the analy- 
sis of three of the four objects. We are aware 
that the results for the four galaxies suffer 
partially from low statistics and incomplete- 
ness due to personal bias and observational 
constraints (such as the beamsize). However, 
we feel that these effects, to first order, affect 
a comparison in a similar way and that it is 
valid to try to find global trends function 
of Hubble type. 

For the sake of brevity, rather than first 
show the statistical properties of the holes in 
IC 2574 separately, we immediately go to a 
comparison with the other three objects. In 



Figs. p!5|-p!9| we compare the distributions of 
swept-up matter, the energies, diameters, ex- 
pansion velocities and ages of the holes in the 
4 galaxies. In each graph we plot in the form 
of histograms the relative number of holes, in 
percent, in order to make a direct compari- 
son possible. In each histogram, the data for 
the individual objects are binned in the same 
way. To improve the presentation, however, 
the respective bins were slightly shifted when 
plotting the results. 



Fig. ^ shows an overlay of the indicative 
masses of the H I holes, binned logarithmi- 
cally. The values range from typically 10^ 
to 5 X 10^ M0. Even though a two-sample 
Kolmogorov-Smirnov test rejects the hypoth- 
esis of pairwise equality of the distributions, 
a casual inspection shows that the indicative 
masses span roughly the same range. The 
mean masses (in lO^M©) for the dwarf galax- 
ies (Ho II: 10.5, IC 2574: 8.4) are slightly 
higher as compared to the spirals (M 33: 4.1, 



20 



M 31: 3.2). This is somewhat surprising if 
one reahzes that the volume density in the 
plane of a spiral galaxy is typically higher by 
a factor of three as compared to that found 
in a dwarf galaxy. H I column densities have 
a much smaller variation from object to ob- 
ject, though (see Tab. ??) which probably ex- 
plains why the swept up matter works out to 
be more or less the same. As a warning, it 
should be born in mind that, in order to de- 
termine the H I densities we smeared out any 
fine scale structure, potentially adding some- 
what to the errors. 
Fig. 
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15, showing in- 



is similar to Fig. 
stead the distribution of the logarithm of the 
energies needed to create the observed H I 
holes (under the assumption of them being 
created by supernova explosions). The de- 
rived energies range from lO^'' to 10^^ ergs. 
We find the important result that the dis- 
tributions are identical, at least to first or- 
der (mean energies (10^^ ergs): Ho II: 1.4, 
IC 2574: 1.1, M 33: 1.6, M 31: 1.5). A two- 
sample Kolmogorov-Smirnov test reveals that 
there is no significant pairwise difference at 
the 95% level between the four energy distri- 
butions. From this we conclude that regard- 
less of Hubble type of the parent galaxy, the 
star clusters (or whatever other process cre- 
ated the H I holes) deposit more or less the 
same amount of energy into the ISM. 

Fig. shows an overlay of the size distri- 
bution of the holes found in the four galaxies. 
In this and the following plots the bins are on 
a linear scale. Note that there is a clear differ- 
ence with Hubble type! The size distribution 
for holes in M 31 and M 33 cuts off sharply 
near 600 pc. In contrast, holes in IC 2574 and 
Ho II reach sizes of 1200 to 1500 pc, respec- 
tively. The lack of holes with sizes smaller 
than ~ 100 pc is due to our resolution limit. 



The fact that the holes get larger for "later" 
Hubble types (i.e., galaxies with lower mass) 
can be understood as follows. Smaller galax- 
ies have lower masses and their mass surface 
density in the disk is lower as well. The one- 
dimensional velocity dispersion varies little, 
if any, among galaxies. We find in the four 
objects under scrutiny, taking into account 
an uncertainty of ± lkms~^, a range of 6-9 
kms~^. This, combined with the lower mass 
surface density, leads to dwarf galaxies hav- 
ing much thicker disks with aspect ratios of 
1:10 rather than 1:100. Hence, for the same 
amount of energy deposited, as argued above, 
an H I shell can grow much larger, because 
of a lower gravitational potential and a lower 
ambient density. Moreover, in dwarf galaxies 
shells are prevented from breaking out of the 
disk because it is thicker (see Sec. |3l^ ). 

The age distribution of the holes (Fig. |T8|) 
also shows a sequence with Hubble type. The 
mean values for the ages of the holes are: 
10 X 10^ years (M 31); 12 x 10^ years (M 33); 
23 X 10^ years (IC 2574) and 62 x 10^ years 
(Ho II). In the spirals no holes are found 
which are older than 30 Myr whereas the dis- 
tribution in dwarfs is flatter and spans an age 
range of up to 120 Myr. This is not unex- 
pected since in more massive galaxies shear 
and mixing due to the passage of density 
waves tend to destroy holes. In dwarfs, un- 
less a shell is hit by a neighboring explosion, 
and subsequently expanding structure, a fea- 
ture can potentially persist for a long time, 
until the rim disperses and mixes with the 
surrounding ISM. 

Finally, the distribution of the expansion 
velocities is presented in Fig. |T^. The pre- 
dominantly low values for the dwarfs (espe- 
cially Ho II) reflect, again, the fact that we 
see many large holes in the dwarfs which have 
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almost stalled. A word of caution should be 
added as the velocity resolution of the IC 2574 
and Ho II observations was a factor of about 
3 better than those of M 31. Very young (i.e., 
smaller than ~ 100 pc) and fast expanding 
holes are missed, of course, but this holds true 
for all galaxies in the sample. 

Figure ^ displays an alternative way of 
looking at the properties of the H I holes 
on a galaxy-wide scale. It shows the expan- 
sion velocities of the shells in the four nearby 
galaxies plotted clS db function of their diam- 
eter (crosses). Based on the hydrodynamical 
models for expanding shells (Chevalier 1974), 
we plotted lines of constant energy for the ex- 
panding structures in Fig. (assuming a con- 
stant density of the ambient medium of 0.15 
cm~^). The plotted curves represent energies 
at a level of 10^^ 10^°, 10^^ and 10^^ g^gg^ 
respectively (from left to right). The evolu- 
tion of an individual hole of a given energy is 
along an equi-energy line from the top of the 
plot downwards. In addition, lines of equal 
ages are shown in the figure, ranging from 10'' 
years (steepest slope) via 5 x 10^ years to 10^ 
years. 

There is a marked difference between the 
two spiral galaxies and the dwarfs. The H I 
shells in M 31 and M 33 are smaller and the 
range of expansion velocities is larger. In fact, 
when looking at the objects in order of de- 
creasing mass (M 31 - M 33 - IC 2574 - Ho II) 
there seems to be a gradually shifting pattern. 
Whereas in the larger spirals there is an ab- 
sence of older holes (older than a few xlO'' 
years, holes in IC 2574 range between 5 x 10^ 
and 10^ years and holes in Ho II cluster near 
10^ years 

This suggests an alternative explanation 
for why we don't see younger holes in IC 2574 
and Ho II. The SF histories of these dwarf 



galaxies could be quite different from the spi- 
rals. In the latter, the global SF rate is, to 
first order, constant with time. Therefore, 
there is always a population of young holes. 
In contrast, in gas-rich dwarf irregulars, SF 
occurs in bursts which transforms them into 
an H II or blue compact dwarf (BCD) galaxy. 
It might be that IC 2574 and Ho II have 
undergone a recent burst, perhaps triggered 
by a close passage with another member of 
the M 81 group, and are now in their post- 
starburst phase and that this is why we don't 
see many young holes. Our data alone don't 
allow us to choose between the two possibili- 
ties. 

Oey & Clarke (1997), using the standard, 
adiabatic shell evolution predict the size dis- 
tribution for the populations of OB superbub- 
bles. In their paper they compare their pre- 
dictions with observations of the SMC and 
the three same galaxies as we compared our 
IC 2574 data with, namely M31, M33, and 
Ho II. Without going into detail, their pre- 
diction fits the observations for the SMC sur- 
prisingly well. This despite the fact that they 
ignore the selection effect (missing short spac- 
ings) which renders the SMC maps insensitive 
to large structures which should dominate the 
appearance of dwarf galaxies. Also, they em- 
phasize stalled shells, i.e., objects which have 
reached their largest diameter. This is in con- 
tradiction with the fact that the SMC shells, 
especially the smaller ones, are all found to be 
rapidly expanding. Their approach works less 
well for the other objects. They ascribe this 
to the data being incomplete, not only when 
considering the smallest detectable shells, but 
overall, due to the reduced sensitivity as com- 
pared to the SMC data. In that respect, our 
IC 2574 are not better than the maps on M31, 
M33 or Ho II. Hence, a thorough test of Oey & 
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Clarke's predictions will have to await higher 
resolution and higher sensitivity data on the 
nearest (dwarf) galaxies. 

6.3. The intrinsic shape of the shells 

One of the results from our analysis is that 
over 90% of the holes are well approximated 
by a circular shape (the few exceptions are 
discussed in Sec. |4.2| ). This is a major differ- 
ence with other galaxies studied thus far in 
similar detail (i.e., Ho II, M 31, M 33, M 101, 
NGC 6946) which tend to have much more el- 
liptical holes. One attempt to understand this 
is in terms of projection effects (different incli- 
nations for the galaxies observed thus far). If 
the shells are intrinsically triaxial, their pro- 
jected shape on the sky will invariably be an 
ellipse with an ellipticity which varies as a 
function of the intrinsic axial ratios, the incli- 
nation of the galaxy, and the position angle of 
the major axis of the hole within the plane of 
the galaxy. This should lead to a wide range 
of ellipticities. The fact that we don't see this 
in IC 2574 can only mean that the shells in 
this galaxy are intrinsically spherical. 

This can be understood by considering the 
shapes of the rotation curves of the galax- 
ies. IC 2574 is the only object which ex- 
hibits a rising rotation curve throughout, as 
is typical for dwarf galaxies. As the rotation 
curve does not turn over, leading to a fiat part 
dominated by differential rotation, there is no 
shear to distort and ultimately destroy any 
large scale structures such as H I shells. For 
a more complete discussion on this topic, see 
Silich et al. (1996) and Palous et al. (1990). 
Indeed, PWBR92 estimate that the amount 
of shear over 1 or 2 kpc in radius, once the 
flat part of the rotation curve is reached in 
Ho II, is sufficient to stretch the holes to an 
axial ratio of 0.5 over a timescale of 10* years. 



IC 2574, however, is the only object studied to 
date which shows solid body rotation almost 
throughout the disk (MCR94) which ensures 
that holes stay circular until long after their 
creation. 

There is one fly in the ointment. Why 
should the shells in IC 2574 be intrinsically 
spherical? The expansion rate perpendicular 
to the disk should be different from that in 
the disk. In the direction perpendicular to the 
disk, the volume density rapidly decreases, fa- 
voring faster expansion in that direction. But 
on the other hand, the gas which moves up- 
wards starts to feel the gravitational pull from 
the mass surface density below it, slowing it 
down. It thus seems as if nature conspires to 
create, at least in this particular object, holes 
which are near spherical in shape. 

Based on the absence of elliptical holes, an 
interesting conclusion can be drawn concern- 
ing magnetic fields. Usually, the influence of 
magnetic fields is neglected in numerical stud- 
ies of expanding bubbles in the ISM. However, 
a few analytical approaches have been put for- 
ward in recent years to investigate the effects 
of the evolution of shell-like structures in the 
presence of magnetic fields (e.g., Tomisaka 
1991). The results of these investigations are 
that if global magnetic fields are present, this 
would inhibit the expansion of a supernova 
driven bubble in directions perpendicular to 
the field lines. The expansion along the field 
lines would not be affected, thus leading to 
elliptical, elongated structures. The fact that 
we observe only spherical holes in IC 2574 
shows that either magnetic fields play a minor 
role or that the magnetic fields in IC 2574 are 
really weak. 
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7. Summary 

High resolution H I and Ha imaging us- 
ing the VLA and the Calar Alto observatory 
of the dwarf galaxy IC 2574 have been pre- 
sented in this paper. The main results and 
conclusions are: 

1. IC 2574 is a gas-rich dwarf galaxy 
which shows a stunning amount of detail in 
the form of H I shells and holes in its in- 
terstellar medium. These features are similar 
to those found in the Galaxy, the Magellanic 
Clouds, M 31, M 33, M 101 and NGC 6946, 
IC 10 and Ho II. From the global H I profile 
we derive a total H I mass of 9.0 x 10^ M0. 

2. Using two independent approaches, the 
thickness of the H I layer of IC 2574 has been 
derived to be ^ 350 pc. This implies that 
IC 2574, as is the case for Ho II, is a 'thick' 
dwarf galaxy (compared to more massive spi- 
rals) which can be understood in terms of its 
lower gravitational potential. 

3. Large scale structures emptied of H I 
(believed to be expanding H I shells) leave 
roughly circular footprints, H I holes. They 
are prominent in the high resolution VLA 
channel maps of IC 2574. In total, 48 H I 
holes have been catalogued which are dis- 
tributed across the entire galaxy. The di- 
ameters of the holes range between 100 pc 
and 1200 pc. The lower limit of the sizes 
is due to the limit of our spatial resolution. 
Most of the holes are found to be expand- 
ing with typical expansion velocities ranging 
between 8 and 12 kms~^. The indicative 
ages of the H I holes range between 10 and 
60 X 10^ years. Typical energies that are re- 
quired to create the observed features are of 
the order of 10-100x10^° ergs assuming that 
they were created by the combined effects of 
stellar winds and supernova explosions. The 



H I masses that were present at the position 
of the holes before creation (and that might 
now be present on the rim of the shell) are of 
the order of 10-100x10^ Mq. 

4. In addition to the well defined holes, 
some large scale (> 1000 pc), coherent fea- 
tures are visible in the channel maps. They 
might be the remainder of an older shell pop- 
ulation. 

5. The radial expansion of the holes, the 
indicative ages and the energy requirements 
for their formation suggest that they have 
been created by the combined effects of stel- 
lar winds and multiple supernova explosions. 
Despite their large dimensions, the energy re- 
quirements are not extreme and infall of ma- 
terial such as clouds is not needed to explain 
the observed features. However, alternative 
mechanisms cannot be ruled out. 

6. The shells can grow to these large di- 
mensions because of several conditions which 
are fulfilled in dwarf galaxies. The volume 
density in the plane is low, which facilitates 
expansion. In the direction perpendicular to 
the disk, the gravitational pull is smaller than 
in a massive spiral. Also, because of the thick 
H I layer, shells are easily contained and un- 
fikely to blow-out. Lastly, solid body rota- 
tion and a lack of spiral density waves prevent 
holes from being rapidly destroyed. 

7. A comparison with other galaxies stud- 
ied in similar detail so far (M 31, M 33 and 
Ho II) shows that the size distribution of 
H I holes found in a galaxy is related to 
its Hubble type in the following way. The 
size of the largest H I shells is inversely pro- 
portional to the global gravitational potential 
(and hence mass surface density). The shells 
in dwarf galaxies show lower expansion veloci- 
ties, most likely due to the fact that the shells 
are preserved much longer due to the absence 
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of shear (differential rotation). The energies 
needed to create these structures, though, are 
found to be roughly the same for all types of 
galaxies. If we follow the conventional view 
that star forming regions are responsible for 
the structures found in H I, the energy output 
of a typical star forming region does not seem 
to be related to the host galaxy. 

8. The fact that virtually all holes in 
IC 2574 are circular in shape is a major differ- 
ence with the other galaxies studied in com- 
parable detail and can be explained, again, 
by the absence of shear and spiral density 
waves in IC 2574. This finding also indicates 
that magnetic fields play only a minor role in 
shaping expanding holes or that the magnetic 
fields in IC 2574 are weak. 

9. Current star formation, as traced by Ha 
emission, is predominantly found along the 
rims of the larger H I-holes, indicating prop- 
agating star formation. A detailed compari- 
son of the Ha and the H I properties yields 
that, on a global scale, star formation is only 
present if the H I surface density reaches val- 
ues larger then 10^^ cm~^ (on a linear scale of 
~ 95 pc). 

In summary, our results on IC 2574 con- 
firm the findings by PWBR92 and show that 
Ho II is not a special case. Our analysis, and 
the comparison with other galaxies studied in 
similar detail suggest that there is a trend 
of hole characteristics with Hubble type (or 
mass) . If the conclusion that the energy out- 
put of a SF region is independent of the host 
galaxy holds then it follows that there is a 
lower mass limit for the formation of dwarf 
galaxies. In other words, if the mass of a 
galaxy falls below a certain threshold, the first 
burst of star formation will tear the whole 
galaxy apart (see also, e.g., Puche & Westp- 
fahl 1994). Mac Low & Ferrara (1998) pre- 



dict that this lower mass limit will be near 
10^ Mq, which is well below the mass of the 
objects dealt with in this paper. 
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Fig. 1.— R-band CCD-image of IC 2574 ob- 
tained with the Calar Alto 2.2-m telescope. 



Fig. 2. — Continuum-subtracted Ha image 
of IC 2574 obtained with the Calar Aho 2.2- 
m telescope. The total exposure time is 45 
minutes. The brightest stars that were not 
fully subtracted by removing the continuum 
were blanked. 

Fig. 3. — Mosaic of every second chan- 
nel map of IC 2574 in which line emission 
was detected. The velocity of each channel 
in kms^^ is plotted in the right hand cor- 
ner. Greyscale range from (white) to 6 
mJy beam^^ (black). Note that the beamsize 
(6.4" X 5.9" as indicated in the box of the top 
left hand panel) is smaller than the resolution 
of this plot. Coordinates are with respect to 
B1950.0 equinox. 

Fig. 4.— Global H I profile for IC 2574 
(squares) obtained by integrating the channel 
maps over the area of the source after sub- 
traction of the continuum background. The 
data have been flux corrected (see text) and 
are corrected for primary beam attenuation. 
The Flux F is given in Jy. The channel spac- 
ing is 2.58 kms^^. For comparison, the pro- 
files published by Rots (1980) (full line) and 
MCR94 (dotted line) are also plotted. The 
arrow indicates the velocity range affected by 
contamination due to galactic emission. 

Fig. 5. — H I surface brightness map of 
IC 2574. The greyscale is a linear representa- 
tion of the H I surface brightness (0 — 2406 K 
km s^^). The beamsize (6.4" x 5.9") is indi- 
cated in the lower left. 

Fig. 6. — Isovelocity contours of IC 2574 of 
the 30" convolved cube. The velocities are 
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given in kms"^. The greyscale image is a lin- 
ear representation of the H I surface bright- 
ness map at full resolution. A circular 30" 
beam is plotted in the lower left hand corner. 



Fig. 7. — Velocity dispersion (second mo- 
ment) map of IC 2574. The greyscale is a 
linear representation of the velocity disper- 
sion ranging between 4 kms^^ (white) and 13 
kms~^ (black). The beamsize (6.4" x 5.9") is 
indicated in the lower left. 

Fig. 8. — Typical examples of position veloc- 
ity diagrams, a) pV diagram of a Type 3 hole 
(hole number 21). b) rV diagram of the same 
hole, c) pV diagram of a Type 1 hole (num- 
ber 30) and d) pV diagram of a Type 2 hole 
(number 33). 

Fig. 9.— Location of the H I holes in IC 2574 
(with numbers) . The greyscale map is a linear 
representation of the H I surface brightness 
map. For the largest holes (13, 23, and 26), 
see text. 

Fig. 10. — Ha image of IC 2574 overlaid with 
the the discovered holes in its neutral inter- 
stellar medium (with numbers). Regions A-C 
are enlarged in Figs. |12|-14, respectively. Note 
that almost all Ha emission in IC 2574 is lo- 
cated on the rims of H I shells (see text). 



Fig. 11. — Ha image of IC 2574 overlaid with 
one contour of the H I surface brightness map 
at 1.7 X 10^^ atoms cm~^. Note that regions 
where we find Ha emission coincide with high 
H I surface brightness. 
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12. — Enlargement of region A from 
The circles indicate where we de- 



Fig. 
Fig. 

tected the H I holes (for identification of the 
H I holes, please refer to Fig. [1^). The boxes 



show how we defined groups of H II regions 
within IC 2574 (with sequence numbers). 

Fig. 13. — Enlargement of region B from 
Fig. The circles indicate where we de- 
tected the H I holes (for identification of the 
H I holes, please refer to Fig. |10|). The boxes 
show how we defined groups of H II regions 
within IC 2574 (with sequence numbers). 

Fig. 14. — Enlargement of region C from 
Fig. The circles indicate where we de- 
tected the H I holes (for identification of the 
H I holes, please refer to Fig. |10|). The boxes 
show how we defined groups of H II regions 
within IC 2574 (with sequence numbers). 

Fig. 15. — Comparison of the relative distri- 
bution, in percentage, of the indicative masses 
of the H I holes in Ho II, IC 2574, M 33 and 
M 31. 

Fig. 16. — Comparison of the relative distri- 
bution, in percentage, of the energies required 
to produce the H I holes in Ho II, IC 2574, 
M 33 and M 31. 

Fig. 17. — Comparison of the relative distri- 
bution, in percentage, of the diameters of the 
H I holes in Ho II, IC 2574, M 33 and M 31. 

Fig. 18. — Comparison of the relative distri- 
bution, in percentage, of the ages of the H I 
holes in Ho II, IC 2574, M 33 and M 31. 

Fig. 19. — Comparison of the relative distri- 
bution, in percentage, of the expansion veloc- 
ities of the H I holes in Ho II, IC 2574, M 33 
and M 31. 

Fig. 20. — The relation between expansion ve- 
locity and diameter in each of the four galax- 
ies. Crosses represent observed values. The 
plotted curves are lines of constant energy 
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(based on Chevalier 1974) at levels of 10^^, 
10^°, 10^^ and 10^^ ergs (from left to right). 
The straight lines correspond to constant age 
and are drawn at 10^ years (steepest slope), 
5 X 10^ years and 10^ years. 
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Table 1: General information on IC 2574 



Object 


IC 2574 


Other names 


UGC 5666, DDO 81, VII Zw 330 


Right ascension (1950.0)^ 


10''24™41!3 


Declination (1950.0)^ 


68° 40' 18" 


Adopted distance^ 


3.2 Mpc 


Inclination'^ 


77±3° 


Scale (1" corresponds to) 


15.5 pc 


Apparent magnitude (B)'^ 


10.8 mag 


Corrected distance modulus 


mB-MB=27.6 


Corrected absolute magnitude (B) 


-16.8 mag 


Blue luminosity Lb 


7.8x10^ Lb0 


HI mass Mhi^ 


9.0x10^ Mq 


HI mass to blue light ratio Mhi/Lb 


1.15 Mq/Lbo 



"Dressel & Condon 1976 

''see discussion in Puche et al. 1992 

'^Martimbeau et al. 1994 

''de Vaucouleurs et al. 1991 (RC3) 

''this study 
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Table 2: Setup of the VLA during the observations 



VLA configuration 
baselines 

date of observation 
total time on source 



B 

1 - 54kA 
(0.21 - 11.4km) 
1992 Jan. 18 
1992 Jan. 19 
647 min 



C 

0.34 - 16 kA 
(0.073 - 3.4 km) 
1992 Apr. 20 

133 min 



D 

0.166 - 4.9 kA 
(0.035 - 1.03 km) 

1991 Mar. 03 

1992 Sep. 03 
245 min 



total bandwidth 
No. of channels 
velocity resolution 
central velocity 
angTilar resolution^ 
linear resolution*^ 
rms noise^ 



1.56 MHz 
128 
2.58 kms-^ 
38 kms~^ 
6.4 X 5.9 arcsec 
100 X 91 pc 
0.7 mJy X beam^-^ 
11 K(^lviu 



■^based on robust data cube; see text 
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Table 3: Observed properties of the H I holes 
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Table 4: Derived properties of the H I holes 



# d 'Vexp age t mass Mhi energy 

(pc) (kms~i) (10^ yr) (10"^ Mq) (lO"'" ergs) 
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Table 5: Cross identification of our grouping of H II regions with the ones defined by MH94 



Number 


"1950.0 


5l9B0.0 


cross identifications 

Miller & Hodge 1994 (MH94) 


1' 


10 


25 20.2 


68 41 40 


246,251 


2' 


10 


25 10.1 


68 40 50 


181,190,193,197,203,206, 










207,218,221,226,227,230 


3' 


10 


25 10.2 


68 40 32 


208 


4' 


10 


25 10.4 


68 40 19 


195,213,214,227 


5' 


10 


25 05.7 


68 41 15 


184 


6' 


10 


25 04.5 


68 39 39 


182,174 


7' 


10 


24 59.5 


68 41 51 


136,143,155,158,165 


8' 


10 


24 54.1 


68 40 53 


125,126,131 


9' 


10 


24 48.9 


68 40 31 


117 


10' 


10 


24 49.9 


68 40 18 


120 


11' 


10 


24 49.5 


68 39 43 


115,121 


12' 


10 


24 48.7 


68 39 27 


113,116 


13' 


10 


24 41.3 


68 38 58 


99,100,106 


14' 


10 


24 44.5 


68 40 07 


111 


15' 


10 


24 42.0 


68 40 37 


105 


16' 


10 


24 40.1 


68 40 05 


94,96,102,104 


17' 


10 


24 .33.1 


68 39 33 


76,80,83,84 


18' 


10 


24 20.4 


68 37 55 


57 58 60-63 65 67 


19' 


10 


24 13.9 


68 37 44 


50 


20' 


10 


24 09.3 


68 37 24 


46 


21' 


10 


23 57.0 


68 38 37 


17,19,20 


22' 


10 


24 02.7 


68 39 12 


30,33 


23' 


10 


24 02.5 


68 39 23 


31 


24' 


10 


24 00.1 


68 39 17 


25 


25' 


10 


24 04.5 


68 39 58 


24,28,32,34,36,38-41,43-46 


26' 


10 


23 59.7 


68 39 58 


24 


27' 


10 


24 49.1 


68 43 30 


118 


28' 


10 


24 58.4 


68 42 29 


144 


29' 


10 


25 05.7 


68 42 35 


186 


30' 


10 


25 03.9 


68 43 04 


178 


31' 


10 


24 59.4 


68 43 23 


148 


32' 


10 


24 56.1 


68 43 24 


133,134,137,138 


33' 


10 


24 58.0 


68 43 52 


140 


34' 


10 


25 03.0 


68 43 51 


167,169,179,180 


35' 


10 


25 07.2 


68 43 26 


192,196 


36' 


10 


25 08.4 


68 43 49 


194,198,202 


37' 


10 


25 12.0 


68 44 06 


210,215,224,225 


38' 


10 


25 14.8 


68 43 51 


235 


39' 


10 


25 17.9 


68 43 51 


237,238,240,242,244,245 


40' 


10 


25 15.7 


68 43 20 


233,237,240 
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Table 6: Derived properties of grouped H II regions 



Number 


F(Ha) 


Phi I 


THI.efF 


Mhi 




PHI 


fHI 


remarks 




[10 erg cm s J 


r^™— 3i 
[cm [ 




[10"" M0j 


km s ^ ^ 




fi n20 — 2l 
[iU cm J 




1 


30.8 


0.7 


201 


14 


16.3 


1.68 


13.9 


Hl-clump 


2 


434.2 


1.8 


225 


24 


26.1 


2.05 


18.9 


Hl-clump 


o/ 


32.9 


0.85 


117 


5.8 


26.3 


3.52 


16.9 


new tmy hole (U<DUpc), rim hole 38 


4 


41.9 


0.8 


148 


10 


19.7 


3.00 


18.3 


diffuse HI 


c:/ 



13.7 


0.5 


152 


7.4 


22.4 


2.05 


12.8 


on rim of hole 39, diffuse HI 


6 


34.5 


0.8 


245 


21 


18.3 


1.39 


14.0 


Hl-clump 


7 


152.4 


0.6 


270 


17 


22.4 


0.84 


9.3 


rla rills hole 32! 


Ot 

a 


55.2 


0.7 


196 


14 


23.7 


1.81 


14.6 


on rim of hole 27 


9 


5.2 


0.8 


111 


4.8 


26.3 


3.41 


15.6 




10 


5.2 


0.6 


107 


4.7 


21.7 


3.73 


16.4 


on rim of hole 30 


11 


40.8 


0.6 


100 


3.9 


28.9 


5.23 


15.6 


TT " riir\r> 

Ha-em. on rim ot hole 28 


1 n/ 
12 


9.9 


0.65 


103 


4.4 


25.0 


3.92 


16.6 


new tmy hole (U<60pcj 


1 o/ 

13 


44.7 


0.6 


295 


28 


23.7 


1.06 


12.9 


" i? 1 1 J'Hr TTT 

on rim ot hole 22, dittuse HI 


14 


8.9 


0.9 


88 


2.7 


26.3 


3.85 


13.9 


on rim ot hole 24, dittuse HI 


1 r/ 

15 


9.9 


0.6 


155 


7.9 


20.0 


2.06 


13.1 


on rim ot hole 24. dittuse HI 


lo 


102.2 


0.9 


163 


10 


34.2 


2.25 


15.0 


on rim ot holes 21,22,24, high Av! 


17' 


86.8 


1.4 


142 


8.5 


25.0 


2.89 


16.9 


on rim of holes 21,22 


18' 


83.8 


0.8 


303 


41 


25.0 


1.43 


18.7 


a lot of HI in region 


19' 


139.0 


1.1 


130 


8.1 


20.4 


3.58 


19.1 


a lot of HI in region 


20' 


13.6 


0.7 


111 


4.3 


22.4 


3.00 


14.0 


a lot of HI in region 


21' 


75.4 


0.85 


233 


17 


18.4 


1.30 


12.5 


new tiny hole (D<60pc) 


22' 


27.5 


1.0 


77 


2.4 


19.6 


5.11 


16.2 


Hl-clump between Ha— em. 


23' 


13.8 


0.8 


95 


3.2 


19.7 


2.37 


14.2 


diffuse HI 


24' 


9.6 


0.7 


93 


2.8 


21.0 


2.43 


12.9 


diffuse HI 


25' 


213.2 


0.8 


270 


23 


31.6 


1.14 


12.6 


Ha same as rim of hole 8! 


26' 


16.2 


0.8 


73 


2.4 


31.3 


6.00 


18.0 


on rim of holes 8,13, high Av 


27' 


35.9 


0.9 


124 


6.2 


25.0 


3.16 


16.1 


on rim of holes 23,29 


28' 


7.8 


0.6 


89 


2.7 


21.0 


3.73 


13.6 


new tiny hole (D<60pc), rim hole 33 


29' 


109.0 


0.6 


93 


4.1 


26.3 


4.96 


18.9 


diffuse HI, on rim of hole 33 


30' 


55.0 


0.8 


136 


7.1 


23.7 


2.75 


15.3 


new tiny hole (D<60pc) 


31' 


120.9 


1.4 


45 


3.3 


34.2 


3.52 


6.5 


hardly any HI, high Av 


32' 


204.3 


1.5 


175 


14 


33.0 


2.54 


18.3 


on rim of holes 31,35, high Av 


33' 


78.8 


0.85 


93 


2.5 


32.0 


3.02 


11.5 


little HI, high Av 


34' 


1004.7 


2.0 


199 


8.9 


25.0 


1.03 


9.0 


little HI, on rim of holes 35,37 


35' 


444.9 


2.0 


154 


11 


23.5 


2.93 


18.5 


a lot HI, on rim of hole 35 


36' 


487.2 


1.9 


183 


17 


28.7 


2.70 


20.3 


a lot HI, rim hole 35, high Av 


37' 


140.4 


1.3 


172 


13 


35.3 


2.48 


17.5 


little HI, rim hole 40, high Av, 


38' 


13.3 


0.8 


98 


5.0 


23.7 


5.17 


20.8 


on rim of hole 42 


39' 


144.8 


1.2 


207 


28 


22.4 


3.07 


26.1 


a lot of HI in region 


40' 


15.4 


0.85 


182 


15 


25.6 


2.42 


18.1 


a lot of HI in region 
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Table 7: Some properties of high mass stars (based on Devereux et al. 1997) 



Type 


Mass (Mq) 


Lbol (L©) 


NLyc 


L(Ha) 


Lboi/L(HQ; 


05 


40 


6.8 X 10^ 


4 X lo^y 


1.4 X 10"^ 


49 


07 


25 


1.0 X 10^ 


4 X 10^« 


1.4 X 10^ 


71 


09 


20 


4.6 X 10'^ 


1.3 X 10^8 


4.6 X 10^ 


100 


BO 


16 


2.5 X lO'' 


2.5 X 10^^ 


89 


282 


Bl 


10 


5.2 X 10^ 


1.9 X 10"^^ 


0.67 


7.8 X 10^ 


B2 


8 


2.9 X 10^ 


4.5 X 10^ 


0.16 


1.8 X 10^ 


B5 


6 


7.9 X 10^ 


4.8 X 10^3 


0.02 


4.0 X 10^ 



Table 8: Summary of H I hole statistics in four nearby galaxies 



Property 


M 3P 


M 33'^ 


IC 2574^ 


Ho 11^ 


Linear resolution (pc) 


100 


55 


95 


65 


Velocity resolution (km s~^) 


8.2 


8.2 


2.6 


2.6 


Average surface density (10^*^ cm~^) 


5 


9 


4 


10 


Average volume density (cm~^) 


0.6 


0.45 


0.15 


0.2 


Average velocity dispersion (km s~^) 


8 


8 


7 


7 


Derived scaleheight (pc) 


120 


100 


350 


625 


Number of holes 


141 


148 


48 


51 


Sensitivity per channel (10^*^ cm~^) 


0.3 


1.0 


0.5 


2.4 



"Brinks & Bajaja 1986 
''Deul & den Hartog 1990 

'^This paper 
''Puclie et al. 1992 
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